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Abstract 



i-^ , It is proposed to use apertures with large acceptance angles to reduce the inte- 

'~^' gration time when studying the emissivity of laser induced plasmas by means of the 
Abel inversion method. The spatial resolution lost due to contributions of angled 

^ . lines of sight to the intensity data collected along the plasma plume diameter is 

i-H I restored by a special numerical data processing. The procedure is meant for the 

^^' laser induced plasma diagnostics and tomography when the integration time needed 

Q ■ to achieve a reasonable signal to noise ratio exceeds a characteristic time scale of 

* ^ ', the plasma state variations which is short especially at early stages of the plasma 

?>^' evolution. 

1 LIBS plasma tomography experiments 

>, 
'nI" . Laser-induced breakdown spectroscopy (LIBS) is widely used for quantitative chemical 

—I. I analysis of an elemental content of various materials [U E]- A typical experimental setup 

>/S I (see, e.g., [31 lU E]) is depicted in Fig. [H A plasma plume is created by a laser ablation 

C^ ' at the distance of two focal lengths from a lens. The lens projects the plasma plume 

image onto the image plane where a detecting device is placed (e.g., a spectrometer 

slit). The plasma plume expands, cools, and radiates. All photons propagating along 

a particular line of sight parallel to a vector e through a point at the distance y from 

k> ; the plasma plume symmetry axis can be collected at the point (0, — ?/,4/), e.g., by a 

Vh ' spectrometer slit. The intensity In{y) measured at the point (0, —y,Af) depends on the 

solid angle Vl in which the light is collected. If fi is small, then the intensity In{y) 

can be approximately viewed as the intensity I{y) in (nearly) parallel lines of sight. In 

this approximation, the measured intensity is related to the plasma emissivity by the 

Abel integral equation that can be solved to recover the emissivity [B]. For this reason, 

spectrometers used in these measurements must have a high /—number in order to restrict 

the solid angle of the acceptance cone. The latter implies, in particular, a small amplitude 

of the measured signal per unit time. In turn, a long integration time is required to 

achieve a reasonable signal to noise ratio. If the plasma plume remained stationary during 



the integration time, the reconstructed emissivity (as well as the plasma temperature 
distribution and local densities of electrons, ions, and atoms) would have corresponded 
to the true (instantaneous) plasma parameters. It appears, however, the plasma plume 
cannot be viewed as stationary during a typical integration time of Ifis or higher. This 
is especially true for early stages of the plasma plume evolution [71 [H] . 

Here the following possibilities are discussed to reduce the integration time. If the 
plasma plume has the spherical symmetry (or, at least it can be approximated as such), 
then it is shown that any device with an arbitrary large acceptance cone can be used to 
collect the data. In particular, spectrometers with low /—numbers can be used. The loss 
of the spatial resolution due to contributions of angled lines of sight into the data function 
In{y) can be eliminated by a proposed numerical data processing such that the intensity 
I{y) in parallel lines of sight is accurately recovered and can be used in the subsequent 
Abel inversion to obtain the plasma emissivity. The data processing is particularly simple 
if the device acceptance cone is slightly restricted by a large aperture of a special shape 
(see Section 3) placed in the lens focal plane (as depicted in Fig. [T]). But with minor 
modifications the reconstruction algorithm is shown to work for any acceptance cone of 
the measuring device (or any aperture). 

If the plasma plume has the axial symmetry, then a measuring device may also have 
an arbitrary large acceptance cone, but the lines of sights that are not perpendicular to 
the plasma plume symmetry axis must be blocked by an aperture in the focal plane whose 
dimension along the symmetry axis is sufficiently small. In other words, only lines of sight 
that lie within a thin slice of the plasma plume (normal to the symmetry axis) contribute 
to the measured intensity. This is similar to the conventional experimental setup (e.g., a 
spectrometer with a high /—number). The difference is that photons propagating along 
all angled lines of sight within the plasma plume slice can be collected to amplify the 
measured signal. A simple numerical data processing algorithm is proposed to accurately 
recover the intensity in parallel lines of sight. 

So, a rational behind the use of larger apertures is a reduction of the integration time. 
This is especially important when studying such a dynamical object as laser-induced 
plasmas. The plasma diagnostics and tomography based on the Abel inversion becomes 
inaccurate when the integration times exceeds a characteristic time scale of the plasma 
state variations which is short especially at early stages of the plasma evolution (e.g., a 
few ns after the ablation process is over). The method proposed below would allow for 
a sufficient signal to noise ratio with integration times in the range of lO's to lOO's ns at 
the price of loosing a spectral resolution. It also allows for the use of spectrometers with 
low /—numbers, which reduces the integration time, but not as much as in the case when 
the spectral resolution is not needed. 
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Figure 1: A typical experimental setup for LIBS plasma measurements. A plasma plume is 
created by a laser ablation at the distance of two focal length from a lens. The light intensity 
is measured in the image plane. An aperture is placed in the lens focal plane. It restricts the 
acceptance cone at each point where the intensity is measured (provided the acceptance cone of 
the light collecting device is wider than that of the aperture) . Photons propagating along any 
particular line of sight parallel to e = e^,,^ that lies within the solid angle of the acceptance cone 
and goes though the point (0,y,0) in the plasma plume are collected at the point (0, — y,4/). 
The measurements give the intensity In{y) as a function of y. For large acceptance cones, In{y) 
cannot be viewed as the intensity in parallel lines of sight and, hence, cannot be directly used 
in the Abel inversion to reconstruct the plasma emissivity. 

2 A generalized Abel equation for apertures with a large acceptance angle 



Suppose that a plasma plume is created on the optical axis of a lens that is at a distance 
of two focal lengths from the plume (see Fig. [1]). Then an image of the plasma plume is 
created in the plane normal to the optical axis that is at the distance of two focal lengths 
but on the other side of the lens. Let the plasma emissivity be given as a function of 
the position vector r, e = e(r). The question to be studied is the amount of plasma 
radiation that can be collected at a point in the image plane. The coordinate system 
is set so that the optical axis coincides with the 2— axis. Consider an infinitesimal area 
element da = dxdz at the point (0, y, 0) (the origin is inside the plasma plume) normal 
to the z— axis. Let e^,^ = (cos yj sin 6', sanpsmO, cos 6) be the unit vector that defines 
the direction of a line of sight through the point (0, y, 0). The geometrical meaning of the 
angles 6 and (p is transparent. The line has the angle 6 with the optical axis. Variations 
of (fi correspond to lines that lie on the cone of the angle 6 whose the axis passes through 
(0,1/, 0) and is parallel to the optical axis. By letting 6 and (f range over the intervals 
[0, 7r/2] and [0, 27r], respectively, all lines of sight through the point (0, y, 0) in the plasma 
plume are parameterized. Without any restriction on the lens diameter, all such lines 
intersect at the point (0, —y, 4/) in the image plane, where / is the lens focal length. 

A lens has a finite size and yet it is subject to optical aberration effects. For this 
reason a diaphragm should be placed in the lens focal plane (see Fig. [1]). The diaphragm 
aperture restricts the "emission" solid angle within which the light can be collected at 



the point (0, —y,4:f). As noted before, spectrometers with high /—numbers are used to 
collect the light (the spectrometer slit is placed in the image plane parallel to the y— 
axis). For a high /—number, the spectrometer acceptance solid angle is smaller than the 
"emission" solid angle. So, not all the plasma radiation that arrives at the observation 
point is collected. If the acceptance solid angle is small, then the measured intensity is 
proportional to the intensity in (nearly) parallel lines of sight going through the plasma 
plume diameter (see below) that is used in the Abel inversion method. 

Let us put aside the question about how exactly the light is collected at the point 
(0, —y, 4/) and assume that all the light coming through the aperture can somehow be 
collected. Alternatively, one can think of any measuring device that has a large acceptance 
cone. As is clear from Fig. [1], any acceptance cone is equivalent to some aperture placed 
in the lens focal plane. So, measuring devices and apertures restricting the "emission" 
cone can be discussed on the same footing. 

The energy collected per unit time, unit frequency, and unit area that was emitted 
into a solid angle dQ in the direction of e^ <^ reads [9] 

Iniv) = / d^cosO / ds£(r0js,y)) = / dQcosO I{9,(f,y) , (1) 

Jn J-oo ^ ^ Jn 

where the dependence on the frequency and time is not explicitly shown (it is not relevant 
for what follows), rQ,^{s,y) = (0,y,0) + e^j^s is the parametric equation of the line of 
sight through the point (0, y, 0) and parallel to the unit vector eg^^, s is the arclength 
parameter, dfl = dOdip sin 9, and the range for the angles 9 and if is determined by the 
shape of the aperture in the focal plane. The function I{9, 99, y) is the light intensity along 
the line of sight parallel to e^^^ through the point (0, y, 0) in the plasma plume. To find 
the dependence of I^ on the aperture geometry, the integration over the solid angle is 
transformed to a double integral over the aperture by a change of variables. Let (m, v) be 
rectangular coordinates in the focal plane such that the origin lies on the lens optical axis 
and the m— axis is parallel to the y— axis (see Fig. [1]). If dA = dudv is the area element at 
the aperture point (m,w), then 

,^ dA cos 9 

^^ = 75^ — ^T~^ 

/^ + M"' + V^ 

and the relation between the angles and rectangular coordinates in the focal plane is given 

by 

V = f tan 9 cos ip , u = f tan 9smip . (2) 

Equation ([T]) can be written in the form 

Uy) = f dA ^ 1(9, ^, y) , (3) 

where the angles 9 and y? in the arguments of / are expressed via u and v by inverting 
relations ^. The integration in ([3]) is carried out over the aperture area denoted by A. 



Consider first the case of a spherically symmetric plasma plume (the case of the axial 
symmetry is discussed in Section 4), i.e., e{r) = e{r) where r = |r|. The distance rQ^^{y) 
between the line r = re(p(s,y) and the origin (the plasma plume center) is given by the 
magnitude of the vector that is the cross product {0,y,0) x ee^^. A simple calculation 
shows that re^^iy) = \ (0, y, 0) x eej = i/(l - sin^ 9 sin^ y?) V2 for y > 0. If J(0, 0, y) = I{y) 
is the intensity along the line of sight parallel to the optical axis, then it follows from the 
assumption of the spherical symmetry of the plasma plume that 

1(9, ip, y) = I{re,M) ' ^e,v{y) = y\/l - sin^^sinV (4) 

because the intensity along any line of sight through the point (0, y, 0) depends only on 
the distance of that line to the plasma plume center. By making use of relations (|2]) to 
express r^,;^ as a function of u and v, Eq. ([3]) can be written in the form 



Uy) = J^ dA (^, ^ f, ^ ^,y l{yKu, V)) , /.(«, V) = ^j /^^,l^, ■ (5) 

If A is an infinitesimal aperture of the area A A centered on the optical axis, then In{y) ~ 
AQI{y) where AQ = AAj p . In this case, the data function can be used as the intensity 
along parallel lines of sight I{y) ~ I^{y)/ AVt that is related to the plasma emissivity e{r) 
by the the Abel integral equation 

It can be solved for e{r). The procedure is known as the Abel inversion [6j: 
1 r^ ^ I'{y) 



roo 

e(r] = -" ' 



roo I (li] poo ^ ^ roo 

- dy r-r^ = / dkkMkr)I{k), I{k) = dye'^yi{y) 

TC Jr Vy ~ f "'0 J-oo 



where I{k) is the Fourier transform of I{y) and Ji, is the Bessel function of index i/ = 0. 
Substituting IQ into the right side of Eq. ([5]), an analog of the Abel integral equation 
(^ is obtained for an aperture with a finite acceptance angle. If Eq. (^ can be solved 
for I{y) for a given Ini^y), then a solution to the generalized Abel equation is found by 
means of the Abel inversion for I{y). In Sections 5 a numerical algorithm is proposed for 
solving Eq. ([5]) in the case of an aperture of a special shape. The case of an aperture of 
a general shape is presented in Section 6. A numerical example of is given in Section 7. 



3 The hyperbolic aperture 

It turns out that there is a particular shape of the aperture for which Eq. (^ is substan- 
tially simplified and becomes convenient for a qualitative analysis of the relation between 
the intensity in parallel lines of sight I{y) and the data function In{y) for large apertures. 





Figure 2: Left panel: A hyperbolic aperture that allows for a factorization in the sum over 
angled lines of sight. It is bounded by the lines v = zbA^ and hyperbolas u = it tan 0^ yjv"^ + /^ . 
Here Ua = ftanOa = /\/l — fJ^'i/fJ-a and u'^ = tan 6a\/ f'^ + A^. For Aq ~ /, the aperture cannot 
be approximated by a rectangle. A rectangular approximation is accurate if A^ <C /. It can 
be used for studies of axially symmetric plasmas as explained in the text. Right panel: An 
illustration to the reconstruction algorithm for a general aperture. Curves 1, 2, and 3 are the 
level curves (hyperbolas) fJ-{u,v) = /i~^, fi{u,v) = ;U^^, and n{u,v) = ^~^_,_]^, respectively, that 
are used to partition the region occupied by the aperture. In each partition region A^f^, the 
function /„(/x) is approximated by a linear function as is explained in the text. Here /Xq is the 
minimal value of the function /x(u, v). 



Let A be symmetric relative to the reflection u — )■ ~u and let A^ be the half of A 
that lies in the half-plane u > 0. Then /^ = 2/^ in the right side of Eq. ([5]). Consider 
the change of variables in the double integral (JSj): {u, v) — )■ (/i, v) where the new variable 
IJ, = fi{u, v) is deflned by the second relation in ([5]). Then dA = dudv = dvdfiJ where the 
Jacobian is J = du/dfi = {dfi/du)^^. The coordinate line fi = const is a hyperbola in 
the {u, t>)— plane: 



u 



fl 



f 



(7) 



It intersects the w— axis at m = fy/T^^Ji^/fi and has slant asymptotes u = \v\y/l — fi'^/fi 
at large \v\. In particular, the line u = (or the w— axis) is the limiting hyperbola as 
H ^ 1. Note that the m— intercept of the hyperbola (JTj) tends to zero, while the slope of 
the slant asymptotes becomes infinite, i.e., the hyperbola tends to the vertical straight 
line through u = 0. Let A^ be bounded by horizontal lines v = ±Aa, the vertical line 
u = 0, and by the hyperbola ([7]) with /x = /Xq < 1 as is shown in the left panel of Fig. [2l 
Then after the change of integration variables, the range of new variables is the rectangle 
—Aft < V < Aa and fJ^a ^ fJ^ ^ ^- A straightforward calculation shows that the double 
integral ([5]) is factorized in the new variables thanks to the "hyperbolic" shape of the 
aperture: 
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dv /■! dfi fi^ 
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Put /i = cos^ so that u = f tan 9a = Ua defines the m— intercept of the hyperbola ([7]) with 
// = /ia that bounds the aperture area A^. Put also taxupa = ^a/ f- Then 

In{y) = 4 sin y^, / " d9 cos^ 9 I{y cos 9) . (9) 

JO 

This equation is simple enough to estimate the signal amplification due to an aperture 
with a large acceptance angle and compare it with signal measurements in reported LIBS 
plasma tomography experiments [5*, "S]. 

In the conventional LIBS plasma tomography experiments, the slit of a spectrometer is 
placed into the image plane parallel to y— axis. However, not all the light that enters into 
the spectrometer slit at the point (0, — y,4/) is collected but rather its only portion that 
illuminates a spectrometer mirror (that further directs the light onto the spectral grating). 
If the aperture in the focal plane is not too small or not present at all, this portion is 
determined only by the /—number of the spectrometer so that the acceptance solid angle 
in Eq. ([1]) is defined by < (^ < 27r and Q < 9 < 9s where tan^^ = l/(2/s) and /^ is 
the spectrometer /—number. Spectrometers in the reported LIBS plasma tomography 
experiments typically have relatively large /—numbers, /^ ~ 5. So, the acceptance angle 
is small, taii9s ^ 9s ^ 10^^. It is then sufficient to calculate the integral (II]) in the 
leading order of 6'^, while neglecting terms of order 9^. In this approximation 7(6', yj, y) = 
I{y) + 0{9'^) (cf. (jl])) and the evaluation of the integral ^ yields: 

Uy) = n9ll{y) + 0{9t). 

The smallness of the factor tt91 in the rate of the energy flux per unit time is exactly 
the reason to have a long integration time to get a sufficient signal to noise ratio, when 
determining I{y) used in the Abel inversion. Suppose now that all the light that comes 
to the point (0, —y,Af) through the hyperbolic aperture is collected. Assuming, for sake 
of simplicity, I{y cos 9) to be a constant /, the evaluation of the integral in ([9]) yields: 

/o^4sin(^,(2^,-sin(2^,))/. 

So the signal amplification (as compared to the spectrometer case) can be estimated by 
a factor of 4sin¥?„(2^„ - sin(2^J)/(7r^2) _ 59 jf q^ ^ iq-i and </?„ ~ ^a ~ vr/4 (i.e., 

Aa ~ Ma ~ /)• 

To increase the angles ipa and 9a-, spectrometers with small /—numbers can be used. 
An effective aperture corresponding to this case is a circular disk of the radius //(2/s). 
If the hyperbolic aperture is fit into this disk, then Eq. ([9]) applies directly, and I{y) is 
recovered by the numerical data processing proposed in Section 5. Otherwise, a circular 
aperture should be treated by the algorithm of Section 6. 

If the spectral resolution is not needed, then a gated ICCD camera can be placed in 
the image plane. The pixel row of the camera is set along the plasma image diameter. The 
light collected by each pixel provides the spatially resolved data In{y)- If the aperture has 



the hyperbolic shape, Eq. (^ apphes directly to recover I{y). Otherwise, the algorithm 
of Section 6 must be used. To study the plasma emissivity in a narrow frequency window 
of interest, a narrow band (or interference) filter can be positioned in front of the detector 
(the pixel array). This method can be particularly useful for diagnostics of laser induced 
plasmas shortly after the ablation process (~ lO's ns). At this stage the plasma exhibits 
a fast dynamics and emits mostly a continuum spectrum radiation. 

4 Axially symmetric plasma plumes 

For axially symmetric plasma plumes, photons propagating along lines of sight that are 
not perpendicular to the plume symmetry axis must be prevented from being collected by 
a detector so that only a plasma plume slice normal to the sjTiimetry axis is visible to a 
measuring device. For spectrometers with high /— numbers, this is achieved automatically 
due to their narrow acceptance cones. But a narrow acceptance cone also eliminates most 
of lines of sight in the visible slice itself that are angled to the optical axis. If i? is a 
plasma plume radius, then the lines of sight through (0, y, 0) in the acceptance cone of 
the angle Og can go through plasma points separated at most by the distance 2Rta,n6s. 
This determines the spatial resolution along the y— axis (as well as along the symmetry 
axis). For Og ~ 10~^, the resolution is such that the data can be collected at about 10 
positions across the plasma plume diameter without a significant overlap of the plasma 
regions from which the light is collected. So, when spectrometers with lower /—numbers 
or pixel rows of an ICCD camera (as suggested above) are used to increase the signal, the 
spatial resolution become even worse along both the directions, the y— and symmetry 
axes. 

The following compromise can be made. The spatial resolution in the direction parallel 
to the symmetry axis is achieved by placing a slit-like aperture extended along the y— axis 
in the focal plane. If a slit of height 2Aa is positioned in the focal plane (— A^ < f < A^ 
in the above notations), then only photons traveling along the lines of sight that are in 
the plasma "slice" of height 2Rtaiiipa = 2RAa/f can reach the image plane. The shape 
of a sufficiently narrow rectangular aperture can well be approximated by the hyperbolic 
aperture with small A^ because the hyperbola ([7]) has the vertical tangent line at the 
intercept point u = Ua (see the caption of Fig. [2]). Thus, Eq. ([8]) or ([9]) remains apphcable 
in this case, but the angle 9a is determined by the measuring device /—number, i.e., 
tan^a = l/(2fs) (the approximation of a small Og is invalid for small fg), and (pa is small 
enough to provide for the spatial resolution along the plasma plume symmetry axis. 

Two-dimensional detectors used in LIBS experiments [U [5] in combination with a 
spectrometer typically have about a few hundred pixels for the spatial resolution along the 
spectrometer slit, while only 10 to 20 spatially resolved data points are taken (each data 
point corresponds to an average intensity over a cluster of pixels) . The algorithm provided 
below allows for an accurate reconstruction of the spatial resolution of just a single pixel 
size. So, it may not be deprived of sense to apply the proposed algorithm to improve the 



spatial resolution even in spectrometers with higher /—numbers. This, however, depends 
very much on the gradients of thermodynamic parameters of the observed plasma. For 
high gradients, this is definitely the case. 

5 The reconstruction algorithm 

Let the light be collected by a linear array of pixel detectors placed along the y— axis in 
the image plane. Each pixel has a size A. The pixel centers are located aX y = yn = nA, 
n = 0, 1,...,A^ (by the symmetry of the observed plasma, it is sufficient to reconstruct 
I{y) (the intensity in parallel lines of sight) only for non-negative y, I{—y) = I{y))- In 
what follows the factor 4 sin ipa in dl]) will be omitted. The energy of all photons collected 
by a pixel n per unit time and unit frequency is given by 

rfl fyn + l/2 J ^ f ^ f^ djji H'^ pn + l/2 j ^ ^ s^ f^ dfi fi 



In = / dylniy)= -1=^ dyIi^^y)= -f^ IM (10) 

"'?/n-l/2 J/la yi-— ft Jyn-1/2 J /la y/ i- — fi 

rf^yn+1/2 
/„(/i) = / dul{u), (11) 

■^mn-l/2 

where yn±i/2 = A.{n ± 1/2) are the boundary points of the nth pixel, and the change of 
the integration variable u = fiy has been carried out. If only the photons propagating 
along parallel lines of sight were registered by the nth pixel, then the data would have 
consisted of values 

/■J/n + l/2 

In= du I{u) = /„(1) . (12) 

Exactly these values must be used in the Abel inversion to reconstruct the emissivity. 
The idea is to find a suitable interpolation for the functions Inilj) using the "sampling" 
points In- Then the integral with respect to /x in fITOl) can be evaluated thus turning Eq. 
(ITO|) into a system of algebraic equations for /„ that can be solved for a given data set /^. 

Let Pn{lj) denote the interval [/i|/„_i/2,/U?/„+i/2] so that PniX) = Pn is the interval 
occupied by the nth pixel. Depending on the values of /i and n, the interval Pn{f^) niay 
have the following positions relative to the intervals occupied by physical pixels: 

PnifJ') overlaps with both Pn-k+i and P„_fc, or PnifJ') is contained in Pn-k (13) 

for fc = 1,2, ...,Nn where A^„ = [(1 — fia)n] + 1 and [x] denotes the largest integer that 
does not exceed x. The integer Nn is defined by the condition yn-N„ < fJ^ayn < yn-Nn+i- 
For a fixed value of n it is not difficult to find the ranges of fi in both the cases flT3l) : 

f^nk< /^ < /^nfc ' or /i-;^+i < ^ < /i+fc , (14) 

l^nk — I 1 /o ' k = 1,2, ..., Nn, /XnNn+1 ~ 1^0- ■ (15) 

For example, for k = 1, the first case in flT3|) holds as long as /^i/„+i/2 remains in Pn, i.e., 
f^yn+i/2 > yn-1+1/2 OT fi > ^ni- ^s yU becomes less than /x^^, the whole interval P„(/i) 



remains in P„_i until the left boundary of Pn{lj) crosses the left boundary of P^-i, i.e., 
l^yn-i/2 > yn-2+1/2 OT fi > fi~2, and SO ou. A Special consideration is required for the case 
k = Nn because, depending on the value of /x^, either both the cases flT3|) . or just the first 
one, or none of them may occur. This issue will be clarified shortly. For time being, the 
values /i^^ are defined by ( TT5l) . i.e., /Xq is assumed not to exceed /U^^+ (both the cases ( TT3|) 
are possible ior k = Nn). 

Consider two integrals J^ dx f{x) and /^, dx f{x) for some f{x) where a < a' < b' < b. 
How can the second integral be approximated by the first one? The following approxima- 
tion is proposed: 

b' - a' [^ 

dxfix)^- / dxfix). (16) 

b — a J a 

The accuracy of this approximation is discussed in Appendix. For a small interval length 
b — a and smooth /, this is a good approximation as is shown in Appendix. If /i is such 
that the interval Pn(/i) is contained in P„_fc, then the approximation (TT6l) applies directly: 

4(/i) = -^ In-k = /i 4-fc , ^nfe+1 < /^ < /^rtfc ■ (l^) 

If fi is such that Pnifi') overlaps with Pn-k and Pn-k+i, then /^(yu) is approximated by a 
linear combination of In~k+i and /„„jt with the coefficients being the lengths of fractions 
of PnifJ') in Pn-k+i and Pn-k, respectively, in units of the pixel size A: 

InifJ') = [l^'O-n- an-k)In~k+l+[K'^- CLn) + an-k)In-k, fJ'nk ^ l^ ^ l^nk , (18) 

where an = n + 1/2. The integral over the interval [fia, 1] in ( fTOl) is split into the sum of 
integrals over the intervals in which the approximations (ITTl) and (ITSl) hold 

The integrals are easily evaluated and the desired system of hnear equations for J„ is 
obtained. There is a simple recurrence algorithm to solve it. Before describing it, the 
issue about the integration limits in the term k = Nn must be addressed. 



Consider the last three integrals in the sum ( fT9|) : 



'""" + /'""" + /'""" , (20) 

in accord with the definition (TTB]) . By examining the possible overlaps of P„(/i) with 
Pn-N„ and Pn-N„+i depending on the value of /i^, the following three possibilities are 
found to occur. If fia > fJ'nN„ where the value of iJ'nNn i^ computed by the rule flTSl) . 
then the last two integrals in f l20|) cannot occur in the sum flTOl) and must be removed. 
If fJ-nNn — k^a < k-nNn where the values of /i^iv^ are calculated by the rule flTSl) . then the 
very last integral in fl20|) cannot be present and must be removed from the sum flT9l) . 

10 



Only if fia < f^nNn all three integrals (12U]) are present in the sum (I2U|) . For algorithm 
implementation purposes, it is convenient to retain the sum (IT9|) as it is. The above three 
cases can be accounted for if, after calculating the values of /i^;, by the rule flTSjl . the 



values of UnNn are redefined by the rule: 

It is straightforward to verify that the unwanted integrals in (1201) (and, hence, in ( IT9ll ) 
always vanish as their upper and lower limits coincide. 

Define the functions 

Fi(a,/3) = £-^l= = l(cos-i(«)-cos-i(/3))+l(av^r^^-/3vT^). 

They are used to evaluate the integrals in the sum (fT9|) . For the central pixel (n = 0), 
one finds 

^0 = -^ — r • (21) 



For n = 1,2, ...,N, Eq. flT9l) yields the recurrence relation 



where 



n n 

Cnl 



Cnk — -Fl(/^„fc; fJ"nk)(^n " -^o(/^„fc5 ^J'nk)^n-k , ^ — 1, 2, ..., iV, 



(22) 



n ) 



-^n = zZ CnkIn-k+1 + 2_. [-^l(/^nfc+l5 /^nfc) ~ ^nk 
k=2 k=l 



n—k 



Note that /„ depends only on Jq, /i,...,/„_i. So by executing the recurrence relation 
with the initial condition (121 p . all the values /„ are found. They can then be used to 
determine the emissivity by a suitable Abel inversion algorithm. There is a vast literature 
on the numerical Abel inversion (see, e.g., [10]), and, for this reason, it will not be discussed 
here. In Section 7, it is demonstrated that for emissivity functions representative for LIBS 
plasmas, the proposed algorithm provides a very accurate reconstruction of /„. Thus, 
the accuracy in the emissivity reconstruction would be determined by the accuracy of a 
particular Abel inversion algorithm. 

6 General apertures 

It is not difficult to extend the proposed algorithm to an aperture of a generic shape. 
Suppose that A is symmetric under the reflection u — ?■ —u and A^ is the half of A for 
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which M > 0. A generahzation to the non-symmetric case is triviaL Integrating Eq. ([S]) 
over the interval occupied by the nth pixel, a generalization of flTOj) is obtained 

^n=2/ ^-4 ^ Uixiu,v)), (23) 

Ja+ (/"' + m^ + f"')^ 

where the function fi{u,v) is defined by the second equation in ([5]). The area A+ can be 
partitioned into regions that lies between the level curves (hyperbolas) of the function 

fl{u, V) = /i+fc , fi{u, v) = /i^^t , 

where ^^j^ are defined in (TTSjl and /!„ = min^^ /i(M,f ). The procedure is shown in the 
right panel of Fig. [2l Let A"^ and A^^ be the regions in which fi^^. < fi{u, v) < fj,~j^ and 
f^nk+i — f^i'^y'^) — l^nki respectively. Then a generalization of Eq. (IT^ to the case of a 
generic aperture reads 






'n = 2E U_ + /.^ ) dA ^^^-^^-^ UKu.v)) . (24) 



The approximations ( TT7|) and (TT8|) are used for /„(/i) in the regions A^^. and A~^, respec- 
tively to evaluate the double integrals in (I21|) . Depending of the shape of A, it can be 
done either analytically or numerically. Put 

/ 2 

F^{B) = 2 f dA ,..\^ ,,, /i"^(M, v), m = 0, 1 , 
for a planar region B. The initial condition (12 ip and the recurrence relation (1221) become 

rfl rH _ jF 

^ - ° . _ ^n ^n n = l,2,...,Ar, (25) 



FM+) ' ai 



where 



C„fc = a„Fi(yl„fc) - a„_feFo(A„fc) , a„ = n + 1/2 , 

-^ri = 2^^'rikln-k+l + 2^\Fl{Ank) — Cnk)In-k ■, 
fc=2 A;=l 

and Anfc is the union of A'^f^ and A~f^. 

7 Numerical results for a hyperbolic aperture 

The reconstruction algorithm has been tested with synthetic data. The emissivity func- 
tion is taken in the form e{r) = eoexp(— r^/(2cr^))(a — cos(6r)) where the values of the 
parameters (arbitrary units) are Eq = 1, a = 3, a = 1.4, and b = 0.5. The emissivity 
profile is typical for the LIBS plasmas observed in experiments ^. The intensity I{y) in 
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Figure 3: Numerical results for a hyperbolic aperture with fi^ = cosOa, Oa = 51.3°. Left panel: 
The solid (black) curve is the graph of the intensity I{y) in parallel lines of sight for the eniissivity 
e{r) described in the text. The dashed (blue) curve is the observed intensity In{v) calculated 
from I{y) as specified in Eq. ([5]). Both the curves are normalized on their maximal values. The 
amplification of the signal is not shown explicitly. It is determined as explained in the text. 
Solid (red) dots are reconstructed values of In from I^. The number of pixels is 60 (A^ = 31). 
Right panel: The same as in the left panel but the number of pixels is 30 (A^ = 16). A small 
inaccuracy of the reconstructed values of In is observed near the maxima of the intensity I{y) 
which occurs through errors in the approximations (J17p and (jlSp . 

parallel lines of sight is calculated by means of Eq. ([6]) . Then the integral intensity over 
all lines of sight within the acceptance solid angle defined by the hyperbolic aperture with 
Ha = COS 9a, 6a = 51.3° is fouud by Eq. ([9]) (the amplification factor isimpa is omitted 
as it is irrelevant for the reconstruction algorithm). The synthetic data are calculated by 
means of Eq. (fTO|) when 60 and 30 pixels are placed across the plasma plume diameter. 
Then the reconstruction algorithm is applied to find the values /„. The results are shown 
in Fig. [HI The left and right panels represent the cases with 60 and 30 pixels, respectively. 
The solid (black) and dashed (blue) curves are the graphs of I{y) and In{y), respectively. 
Both the curves are normalized to their maximal values so that the amplification factor 
due to the vertical size A^ of the hyperbolic aperture cannot be seen. This normalization 
would not make sense for a generic aperture because there is no factorization in the sum 
over angled lines of sight for apertures of a general shape (see Section 6). Yet, the factor- 
ization due to the hyperbolic aperture leads to a mild deviation of In{y) from I{y) when 
they are normalized on their maximal values. The solid (red) dots show the reconstructed 
values of J„. The reconstruction is very accurate for the case of 60 pixels. When the 
pixel size increases, the approximations ( !T7|) and ( !T8|l becomes less accurate as proved in 
Appendix. A reconstruction error increases near the maxima of I{y) when the pixel size 
is doubled (the right panel of Fig. [3]). For a typical pixel detector, the pixel size is about 
0.02 mm. So, for a typical LIBS plasma plume diameter of 2 m,m, there are about 100 
pixel collecting the data. 
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A. Appendix 

Here the accuracy of the approximation flT6l) is assessed. For a continuous function f{x) 
there exist points a < x^, < 6 and a' < x'^ < b' such that 

' dxfix) = {h- a) fix,) , I'= [' dxfix) = {h' - a') fix:) . 

■J a' 

It follows from these relations that 

i' = ^i+A'{f{x:)-f{x,)), 

where A = 6 — a and A' = b' — a'. By the mean value theorem, the last term in this 
equality can be transformed so that 

/' = f / + A'(x:-x.)/'(0 = ^/ + 7 
for some ^ between x'^ and x*. Therefore the error 7 is bounded by 

|7|<MA'A, M = max|/'(x)|. 

[a,b] 

If A ^^ 0, while A'/ A remains finite, the error decreases as A^. 
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